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ABSTRACT 

The theory of gamma ray line emission from solar flares is 

reviewed and revised. It is shown that the recently observed (Chupp 

et al., 1972) line emissions at 0.5, 2.2, 4.4 and 6.1 MeV are due to 

positron annihilation, deuterium de-excitation following neutron 

capture on hydrogen, and the deexcitation of excited states in carbon 

and oxygen. From the observed relative line intensities it is possible 

to determine the spectrum of accelerated protons in the flare region. 

This spectrum is found to be very similar to that of charged particles 

from the flare observed near earth. The total numbers of protons 

at the sun is deduced from the observed absolute line intensities for 

various interaction models. It is found that if the protons at the 

sun have a spectrum which is an exponential in rigidity, the total 

28 

energy in protons is a few times 10 ergs if the gamma rays are 

30 

produced by protons moving down into the sun; and about 10 ergs if 
the gamma rays are produced at the site of the acceleration. 


I 
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INTRODUCTION 

Solar gamma-ray line emissions at 0.5, 2.2, 4.4 and 6.1 MeV 
were detected during the flare of August 4, 1972 by a gamma-ray 
monitor flown on 0S0-7 (Chupp et al. 1972). Line emissions at 0.5 and 
2.2 MeV were also detected on August 7, 1972, but only upper limits 
could be set on the 4.4 and 6.1 MeV lines from this flare. In previous 
papers (Lingenfelter and Ramaty 1967, paper I; Cheng 1972), the theory 
of nuclear reactions in solar flares was treated in detail and pre- 
dictions were made as to the expected fluxes of gamma rays and high 
energy neutrons at earth from such reactions at the sun. Following 
the discovery by Chupp et al. we have reviewed and updated these 
calculations including more recent nuclear cross sections. By 
comparing the predicted emissions with the observations, we can show 
that the observed lines at 0.5, 2.2, 4.4 and 6.1 MeV are produced, 
respectively, by positron annihilation, deuterium deexcitation follow- 
ing neutron capture on hydrogen, and the deexcitation of the first 

12 16 

nuclear levels of C and 0 . Furthermore, from the comparison of 

the calculated and observed line intensities we can deduce the spec- 
trum of the accelerated particles at the sun independent of the assumed 
interaction model. The total number of accelerated particles 
required to produce the observed line emission, however, does depend 
on this model. In the subsequent treatment we shall use two limiting 
models. A thick-target model in which the accelerated particles move 

A ' ' 

from the flare region downwards into the Sun undergoing nuclear inter- 
actions as they slow down in the solar atmosphere; and a thin target 
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model in which the spectrum of accelerated particles is not modified 
during the time in which the nuclear interactions take place. The 
latter model assumes that either the total path length traversed by 
the particles at the sun is small in comparison with their interaction 
length, or that the particle energy loss from ionization and nuclear 
interactions is just balanced by energy gains from acceleration. 

We have recalculated the fluxes of various gamma ray lines 
expected from accelerated particle interactions in solar flares: 
at 0.511 MeV from positron annihilation, at 2.23 MeV from neutron 
capture on hydrogen and at 1.63, 1.99, 2.31, 4.43, 5.5 and 6.14 MeV 
from deexcitation of nuclear levels in C , 0, N and Ne. These calcu- 
lations are based on an ambient solar composition given by H:He:C:N:0 
1: 10"^: 5. 3xl0 -< S 10"^: 9. 2x10"^ (Cameron 1967), and an accelerated 
particle population consisting of protons with spectrum 

N(P) = P 0 _1 exp(-P/P G ). (1) 

Here P is rigidity, P Q is a characteristic rigidity which we treat as 
a free parameter, and N(P) is the total number of protons per unit 
rigidity. In the calculations N(P) is normalized at the sun to 1 
particle of rigidity greater than zero. 

For the thick-target model, the yield of secondaries from a 
particular type of interactions is given by 

Q s = T| f dP N (P) P dx' ct(P') exp (- (x-x ' ) /L) , (2) 

- o ,J o 

where 7] is the number of target atoms per gram of solar material, 
a is the cross section as a function of rigidity, and x and L are the 
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stopping range and nuclear interaction length of protons of rigidity 
P (both measured in g cm"^) . This yield has the units of secondaries 
per incident proton of rigidity greater than zero. 

For the thin-target model, the instantaneous production rate of 
secondaries from the same type of interaction is 

q = n c f dP N (P) 0 CT(P) ( 3 > 

b ■'J o 

-3 

where n is the number density of the ambient solar material in cm , 
and c 3 is the velocity of a proton of rigidity P. The units of q s 
are secondaries per second per proton of rigidity greater than zero. 
PHOTON PRODUCTION 

Line emission at 0.511 MeV is produced from the annihilation 
of positrons. The principal source of positrons in solar flares are 
nuclear reactions of accelerated particles with the ambient solar 
atmosphere. These reactions produce tt + mesons and a variety of 
radioactive isotopes which decay by positron emission. The main 
positron emitters, their formation reactions, threshold energies, 
half-lifes and maximum positron energies are given in Table 1. Except 
for the reactions N 14 (p,n) 0 14 (e.g. Andouze, et al. 1967) and 
N 14 (p,d) c" (Jacobs et al. 1972), the cross sections for these 
reactions were given in paper I. The resultant positron yields are 
given in Figures 1 and 2, for the thick- and thin-target models 
respectively. 

The intensity of the 0.511 MeV line depends on the number of 
positrons that annihilate at the sun. For the thick-target model it 
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is reasonable to assume that all the flare produced positrons will 

ultimately annihilate at the sun. For the thin-target model, however, 

it is possible that a significant fraction of the positrons will escape 

from the sun before they annihilate. The determination of this 

fraction, however, is beyond the scope of the present paper and we 

defer the discussion on this for future research. For simplicity, in 

the subsequent calculations we assume that all flare produced positrons 

annihilate at the sun, but we keep in mind that this assumption may 

lead to too many 0.511 photons in the case of the thin-target model. 

The annihilation radiation yield at 0.511 MeV also depends on the 

mode of positron annihilation. Positrons annihilate either directly 

with a free electron or in a bound state of positronium. In the 

latter case, the annihilations proceed from a state leading 

3 

to two 0.511 MeV photons, or from a S state leading to a 3-photon 

continuum. Since the 3-photon continuum does not contribute to line 

3 

emission at 0.511 keV, and because the probability of forming the S 
state is 3 times that of the state, positronium formation and its 
annihilation produces on the average only one 0.511 photon per positron 
as compared to two photons for free annihilation. This point was 
discussed in detail by Stecker (1969) for positron annihilation in 
interstellar space. However, in the much higher density and temperature 
of the solar atmosphere, we do not expect that positrons will annihilate 

3 

from the S state of positronium, because the collision frequency of 
positronium with ambient electrons is sufficiently high to either 
dissociate the positronium or to cause a transition to the state 
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before the annihilation of the S state. By allowing 2 photons per 
positron at the sun, the time-integrated 0.511 MeV line intensity at 
earth for the thick-target model is 

F(0.511) = -i-2 [Q s (tt + ) + Q s (C " ) + Q s (N 13 ) + Q s (0 14 ) + Q s (0 15 )] (4) 

2rrR 

where the various Q s 's are given in Figure 1. The units of F are 
_2 

photons cm per proton of rigidity greater than zero. For the 

thin-target model we obtain an equation similar to equation (4), with 

F and Q s replaced by the instantaneous rates 0 and q g . The units of 
-2 -1 

0 are photons cm sec per proton of rigidity greater than zero. 

The resultant 0.511 MeV photon fluxes are given in Figures 3 and 4, 
for the thick and thin-target models respectively. 

Next we consider the 2.23 MeV line resulting from the deexcitation 
of deuterium following neutron capture on hydrogen. The principal 
neutron producing reactions and their cross sections were discussed 
in paper I and references therein. To these reactions we have added 
the reaction N^ 4 (p, n) O'*' 4 mentioned above which is important for low 
energy protons (< 15 MeV). The neutron producing reactions are 
summarized in Table 2. The total neutron yields for the thick- and 
thin-target models are then obtained by using the appropriate cross 
sections in equations (2) and (3). Having obtained the neutron yields, 
we have to consider the propagation of neutrons in the solar atmosphere, 
since 2.23 MeV photons are not produced if the neutrons escape from 
the sun or decay in the solar atmosphere. 

Neutron propagation in the solar atmosphere is determined 

predominantly by scattering (elastic and inelastic) and capture by 
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protons. The total neutron-proton scattering cross section a s is 
essentially constant at 20 barns for neutron energies from 1 eV to 
10^ eV, and then drops to about 0.1 barns at 100 "MeV (Hughes and 
Schwartz 1958). The capture cross section, is inversely proportional 
to velocity and is given by ct c — 2.2x10 8~ , where a c is measured 

in barns and c6 is the velocity of the neutron. Since for all energies 
of interest <j a » a , the neutron mean free path in the solar atmos- 
phere is determined principally by neutron-proton scattering. Using 
the columnar density of the solar atmosphere as given by Allen (1963), 
we find that a 10 MeV neutron moving radially outward from the sun 
will probably escape from the solar atmosphere before making even 

17 -3 

one collision if it is produced in a region of density less than 10 cm 
This density corresponds to a depth of about 300 km below the base of 
the chromosphere. Since solar flares occur in the chromosphere or 
corona, it is reasonable to assume that for both interaction models, 
all upward moving neutrons are going to escape from the sun. The 
downward moving neutrons on the other hand decay or are captured, 
depending on whether the capture time is greater or smaller than the 
half-life of the neutrons. From the capture cross section cr c5 the 
capture time t c is given by 


t 


c 


-1 

(c8na c ) 


19 

1.5x10' _ 
n(cm”^) 


sec 


(5) 


As can be seen, t is independent of energy and is inversely proportional 
to the density n of the region where the neutrons interact with the 
ambient medium. From the discussion above it follows that downward 
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directed neutrons probably will not collide with the ambient gas till 

they reach the photosphere where n is of the order 10 cm . Thus 

t c =“ 150 seconds and since the half life of the neutron is 720 seconds, 

it follows that most of the neutrons will be captured before they 

decay. Furthermore, since the neutron capture occurs at a columnar 
24-2 , -2 

depth of about 10 cm or 1.6 g cm while the stopping range of a 

-2 

2.2 MeV photon is about 25 g cm , all Ihe upward moving photons 
resulting from deuterium deexcitation will escape from the sun. 

Assuming isotropic production of neutrons in the interaction region, 
the 2.23 MeV line intensity at earth, for the thick- and thin-target 
model respectively are 

F(2. 23) = ~ Q s (n) (6) 

0(2.23) = ^ 2 q S (n) ’ (7) 

where Q s (n) and q s (n) are the total neutron yields as calculated from 
equations (2) and (3). 

The calculations of the other intensities at 1.63, 1.99, 2.31, 
4.43, 5.5 and 6.14 MeV in Figures 3 and 4 is straightforward. Unlike 
the 0.511 and 2.23 MeV lines, these emissions are prompt, i.e. the 
excited states or secondary products decay by photon emission in a 
time scale much shorter than any characteristic time of the flare 
process . 

Line emission at 4.43 and 6.14 MeV results from the deexcitation 

1 9 i f. 

of the first nuclear levels of and 0 1 , respectively. The inten- 
sity of the 6.14 MeV line is the same as calculated in paper I. The 
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intensity of the 4.43 MeV line is about 50%, greater than in paper I, 
because it is possible to produce q^^( 4.43) the spallation of 0^, 
a process which was neglected in paper I. The calculations for this 
process given in Figures 3 and 4 are based on cross sections measured 
by Zobel et al. (1968). 

Radiation at 2.31 MeV corresponds to the first nuclear level of 

N^. The 2.31 MeV line is produced by the direct excitation of the 

first and second levels of . The latter is at 3.94 MeV and it 

deexcites 96% of the time through the first level, thereby producing 

a 1.63 MeV photon in addition to the 2.31 MeV photon. In addition to 

direct excitation, the first level of can also be populated by 

the decay of 0^ which is produced by the reaction N^(p,n)0^. 

0 1Zh beta decays 99.4% of the time to ^14(2. 31) anc j hence each beta decay 

is accompanied by a 2.31 MeV photon. As mentioned above, the 1.63 

14(3.94) 

MeV line is formed by the deexcitation of N " .In addition, this 

20 

line is also produced by the deexcitation of the first level of Ne 

Finally the 1.99 MeV line results from the deexcitation of C* which 

12 

is formed from the spallation of C . The cross section for this 

process is given by Zobel et al. (1968). 

The principal reason for showing the 1.99 MeV and the 2.31 lines 

is that for detectors of poor energy resolution these lines could in 

principle be confused with the deuterium deexcitation line. As can 

ii a. 99) 

be seen from Figures 3 and 4, however, the C v ' MeV line is 
negligible in comparison with the deuterium line. On the other hand, 
the 2.31 MeV line from could compete with the 2.23 MeV line at 



low P Q ' s . If all the emission at ~ 2.2 MeV is from N^, however, the 
intensities of the 4.43 MeV and 6.14 MeV lines should be about the same 
as that of the 2.2 MeV line. As will be discussed below, this was not 
the case for the August 4, 1972 event. Therefore, for this flare at 
least, we conclude that the 2.2 MeV line is produced almost entirely 
by neutron capture. 

We have summarized in Table 3 the principal mechanism leading to 
line emission in the solar atmosphere. In addition to the lines given 
in Figures 3 and 4, we have listed in Table 3 lines at 5.2 MeV from 
0 15 and and at 7.12 MeV from 0^. From the cross sections of 
Zobel et al. (1968), the intensities of both these lines should be 
approximately 50% of the 6.14 MeV line intensity. 

DISCUSSIONS 

Let us now compare the results of our calculations with the 
observations of Chupp et al. (1972) for the August 4, 1972 flare. 

We defer the discussions on the August 7, 1972 observations, since 
the gamma rays from this flare were observed only after the flare 
maximum and hence they require a more detailed treatment of the 
time dependence of the gamma ray intensities. 

In Figures 5 and 6, the shaded areas represent the time averaged 
ratios of the measured 0.5, 4.4 and 6.1 MeV line intensities to the 
measured 2.2 MeV line. The curves represent the calculated ratios 
as functions of the characteristic rigidity P Q , for the thick-target 
model in Figure 5 and the thin-target model in Figure 6. As can be 
seen, the calculated 4.43 and 6.14 curves are strong functions of P Q . 
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Therefore, the comparison of these ratios with the measurements allows 

us to deduce the value of the characteristic rigidity P Q . We find 

that for both models P 0 has to be in the range 70 to 80 MV. This range 

of P D 's should be compared with values of P Q as obtained from charged 

particle observations near earth. 

According to Bostrom et al. (1972) the peak proton intensities 

6 -2 -I 

after the flare of August 4, 1972 were: j (> 10 MeV) = 10 cm sec , 
j(> 30 MeV) = 2.6xl0 5 cm“ 2 Sec _1 , and j (> 60 MeV) = 8xl0 4 cm" 2 sec _1 . 

From these intensities we can deduce the local proton density u in the 
■ interplanetary medium. If u(P) «: .exp (- P/P 0 ) and if the protons are 
non-re la tivis tic 

u(> P) = j(> P) (mc/e) (P+Pq) -1 . ( 8 > 

From equation (8) and the integral intensities j (> P) given above, 
we can calculate values of u(> P) for various P D 's. For P D = 65 MV 
we obtain proton densities as follows: 

u (> 10 MeV) = 1.53xl0“ 4 , u(> 30 MeV) = 2.67xlp" 5 , and u(> 60 MeV) = 
6.13xl0~ 6 . These numbers are plotted in Figure 7 as functions of P. 
The straight lines are exponentials in rigidity with P Q = 60 MV and 
P Q = 70 MV and they bracket the observed proton densities. The 
charged particle observations at the peak of the proton event are 
therefore consistent with a characteristic rigidity of about 60 to 
70 MV. Furthermore, this characteristic rigidity remains essentially 
the same for the remainder of the particle event (J. King, private 


communication) . 
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The similarity between the proton spectrum as observed near earth 
and the proton spectrum at the sun as deduced from the gamma ray 
observations, seems to imply that, at least for the August 4, 1972 
flare, both the escape of particles from the flare region and their 
propagation in the interplanetary medium are essentially independent 
of energy. 

Let us now calculate the total number of protons at the sun. In 
the case of the thick-target model this calculation is independent 
of the ambient density n but requires the knowledge of the time-inte- 
grated photon flux from the flare. According to Chupp et al. (1972) 
the average intensities of the 2.2 MeV line from 0626 UT to 0632 UT 
on August 4, 1972 was 0.22 photons cm" 2 sec -1 and thus the total photon 
flux in the 6 minute time interval was 80 photons cm" 2 . This is only 
a lower limit because 0S0-7 went into earth eclipse at 0632 before the 
termination of the gamma ray event. The total gamma ray flux, however, 
was probably not much larger, since as indicated by the microwave 
data (Toyokawa observatory, private communication), the acceleration 
of the charged particles probably ceased at 0635 UT, 3 minutes after 
the eclipse of 0S0-7. 

_o 

For a flux of 80 photons cm and P Q ranging from 60 MV to 70 MV, 
we find from Figure 3 that N(> P) is between 1.4 and 6xl0 33 . These 
numbers should be compared with the number of protons released from 
the flare based on measurements of the proton flux near earth. This 
number can be obtained from the local densities shown in Figure 7 if 
we assume that the particles fill some volume V in the interplanetary 
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medium to that density. A conservative estimate of this volume would 
be to assume that it is a cone of opening angle 30° with vertex at 
the sun and height 1.5 A.U. This is consistent with the fact that 
the August 4 flare was located at 15°N latitude on the sun and that 
some charged particles were observed on the Pioneer 10 space probe 

at 2 A.U. from the sun (B. J. Teegarden, private communication). The 

39 3 35 

volume V is therefore ~ 10 cm and from Figure 7, N(> P) is 7x10 

and 1.5x10^ for P Q equal 70 MV and 60 MV respectively. When compared 
with the values of N(P > 0) deduced above from the gamma ray obser- 
vations for the thick-target model, we see that possibly not more 
than about 17„ of the flare accelerated protons could have escaped 
downward into the sun. However, we should note that this conclusion 
is stron^Ly dependent on our estimate of the total number of protons 
at the sun based on the local proton observations. Clearly a more 
detailed understanding of the source function of flare protons in the 
interplanetary is required. 

In case of the thin-target model, the instantaneous photon flux 

at earth directly determines the emission measure at the sun, i.e. 

the product nN(P > 0). From Figure 4 we find that the observed flux 

of 0.22 photons cm ^sec ^ implies that nN is 1 to 2 x lO^cm ^ for P G 

ranging from 70 MV to 60 MV. Then using the values of N(P > 0) as 

deduced from the proton flux near earth, we get that n is 1.5 to 3 x 
9 -3 

10 cm . As before, we note that this density is subject to the 
uncertainty in the total number of protons at the sun obtained from 
the local observations. 



13 


Let us finally evaluate the total energy in the flare accelerated 

protons that produce gamma rays at the sun. For the spectrum given 

in equation (1) and nonrelativistic protons, the average proton 
2 2 

energy is P Q /Me . Therefore, for the thick-target model, the energy 
in the 1.4 to 6x10^ protons as deduced above is 1.2 to 3.7 x lO^ergs. 

For the thin-target model, both the instantaneous gamma ray 
production rate and the instantaneous proton energy loss rate depend 
on the product of the ambient density n and the total number of protons 
at the sun. The ratio of these two rates, therefore, is independent 
of both n and N and depends only on the spectrum of the accelerated 
particles, i.e. P Q in our calculations. This point was first made 
by Lingenfelter (1969). For the proton distribution in equation (1) 
normalized to 1 proton of rigidity greater than zero, the instan- 
taneous energy loss rate W to ionization, excitation and nuclear inter- 
actions in ambient hydrogen of unit density is given in Table 4. As 
can be seen, for P 0 's between 60 to 70 MV, W is very closely equal 
to 1.7x10 °ergs sec . Using the values of Nn given above, we see 
that gamma ray production by flare accelerated protons is accompanied 
by the dissipation of 1.7 to 3.4 x 10 27 ergs sec -1 . For the time interval 

of 6 minutes during which gamma rays were observed, the total dissipated 

30 

energy is 0.6 to 1.2 x 10 ergs. 

30 

This energy is comparable to the total of ~ 2x10 ergs emitted 
in H-alpha by the flare of August 4, 1972 (H. Zirin, private communi- 
cation) and is consistent with the suggestion of Gordon (1954) that the 
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optical energy from flares results from ionization losses of accelerated 
particles. On the other hand, in the thick-target model, the 
ionization losses of the gamma-ray producing protons is only a few 
percent of the observed optical energy; in this case additional energy 
loss may be expected from electrons. 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


- Yield of positron emitters at the sun for the thick- 
target model. 

- Production rate of positron emitters at the sun for 
the thin-target model. 

- Time integrated photon fluxes at earth in the 
thick-target model. 

- Instantaneous photon fluxes at earth in the 
thin-target model. 

- Relative line intensities for the thick-target 
model 

- Relative line intensities for the thin-target 
model . 

- Proton densities in the interplanetary medium 


from local particle measurements. 



Positron Emitters 

Emitter and Decay Mode I Maximum Positron I Half-life Production Threshold Energy 
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TABLE 2 

Neutron Producing Reactions 


Reaction 


Threshold Energy 
Mev/Nucleon 


H^p 

, nn + ) H 

292.3 

He 4 

3 

(p,pn) He 

25.9 

He 4 

(p,2pn) H 2 

32.8 

He 4 

(p,2p2n) H 1 

35.6 

c 12 

( P i n . . . 

19.8 

N 14 

(p,n. . . 

6.3 

o 16 

(p,pn. . . 

16.5 



(p,pn. . . 


17.7 


Gamma Ray Line Emission Mechanisms 
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TABLE 4 

Energy loss rate of the proton distribution in 
Eq. (1) in solar material of n = 1 cm -3 


P (MV) 
o 

20 

30 

40 

60 

80 

100 

120 

200 


W(erg sec “*") 


2 x 10 


-18 


2 x 10 


-18 


1.9 x 10 


-18 


1.8 x 10 
1.6 x 10 
1.5 x 10 
1.4 x 10 
1.2 x 10 
1.1 x 10 


-18 

-18 

-18 

-18 

-18 

-18 


300 
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Figure 5 
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Figure 6 
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Figure 7 



